Spin models on Platonic solids and asymptotic freedom 
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We consider a two-dimensional a-model with discrete icosahedral/dodecahedral symmetry. We present high- 
precision finite-size numerical results that show that the continuum limit of this model is different from the 
continuum limit of the rotationally invariant 0(3) tj-model. 



Recently, there has been interest in the critical 
behavior of two-dimensional cr-models in which 
the spins take values in some discrete subset of the 
sphere. In particular, two groups studied 
the nearest-neighbor a-model 

if = ^ <T, • (1) 

<l]> 

in which the spins have unit length and belong to 
a Platonic solid, i.e. to a tetrahedron, cube, oc- 
tahedron, icosahedron, or dodecahedron. Several 
quantities have been computed, the renormalized 
two-point function, the current-current correla- 
tion function, the finite-size scahng (FSS) curve 
for the second-moment correlation length, and 
the four-point renormalized coupling. Surpris- 
ingly enough, the results for the icosahedral and 
the dodecahedral model are very close to the 0(3) 
ones, suggesting that these three models might 
have the same continuum limit. Patrascioiu and 
Seller [|l]-^ considered these results as evidence 
for the 0(3) cr-model not being asymptotically 
free, since the discrete-symmetry models have a 
finite P phase transition, which cannot be de- 
scribed in perturbation theory. However, the 
overwhelming evidence we have collected in the 
years in favor of asymptotic freedom made Hasen- 
fratz and Niedermayer suggest that, may be, 
the icosahedral and the dodecahedral models have 
an asymptotically- free continuum limit, in spite 



of the fact that the critical point is at a finite 
value of /3. 

In Rcf. 1^ we showed that, by using 
some standard assumptions, the pcrturbative 
renormalization-group (RG) approach predicts 
that the suggestion of Hasenfratz and Nieder- 
mayer cannot be true. If the continuum limit 
of the 0(3) cr-model is correctly described by 
the pcrturbative RG, then any discrete-symmetry 
model cannot belong to the same universality 
class of the 0(3) cr-model. 

The argument goes as follows. Consider the 
Hamiltonian 

H = -f3j2 • '^J -hY,Ini'^^). (2) 

<i,j> i 

where cr^ is an 0(3) unit spin and In{cri) is a poly- 
nomial in (Ti with the following properties: it has 
0(3) spin n; the maxima (or minima) of /„(crj) 
correspond to a Platonic solid; it is invariant un- 
der the discrete-symmetry group of the solid. For 
all Platonic solids, it can be shown explicitly that 
such a polynomial exists. The model (||) interpo- 
lates between the 0(3) {h = 0) and the discrete- 
symmetry model {\h\ — +oo). Now, with quite 
standard assumptions, one can show that Inicr) 
is a relevant perturbation of the 0(3) fixed point. 
In other words, any arbitrarily small perturba- 
tion with discrete symmetry of the 0(3) cr-model 
drives the system to a different fixed point M . 
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Figure 1. FSS function for the spin-1 susceptibil- Figure 2. FSS function for the spin-1 susceptibil- 

ity xi- ity Xi- Here, we restrict the horizontal range to 

1 < a; < 2. 



This argument together with the numerical re- 
sults of Refs. puts asymptotic freedom on 
a dangerous ground since it shows that the con- 
ventional scenario is wrong if the icosahedral or 
the dodecahedral model have the same continuum 
limit of the 0(3) model. We have thus decided 
to extend the previous numerical work and in- 
deed, we have found good evidence that the 0(3) 
model and the icosahedral model belong to dif- 
ferent universality classes: the conventional sce- 
nario is saved. However, the surprising fact is 
that these differences appear only very near to 
the critical point, i.e. for ^ 10^! 

In the numerical simulation we have considered 
the Hamiltonian (||) with 



a', -bat (al+al)+5a', fa^ + a^' 



+2a,a, {at - lOala^ + a^) , (3) 



and h — 0.1. Such a polynomial is invariant under 
the rotation group of the icosahedron and of the 
dodecahedron. We measured the second- moment 
correlation length as defined in Refs. [@,D, and 
the spin-n susceptibilities 



Xn 



^ (Pn (cto • cr^)) , 



(4) 



where Pnix) is a Legendre polynomial, for n = 
1,3,4. Of course, xi is the standard suscepti- 
bility. For each observable 0{L,l3), we consid- 
ered the so-called step function, i.e. the ratio 
0{2L, i3)/0{L, /3), which, in the continuum limit 
should become a universal function of ^(L, f3)/L, 



OiL,P) L 



+ 0(L-,r")- (5) 



We measured the step function of these observ- 
ables in the icosahedral theory (i.e. keeping 
h = 0.1 fixed) and in the 0(3) model, thereby 
extending the results of Refs. |^,^. If the two 
models have the same continuum limit, the func- 
tion computed for h = 0.1 and h = should 
coincide. 

In Figs. |l|, ||we report the numerical results for 
the spin-one susceptibility. The continuous line 
is a fit to the 0(3) data, while the points refer 
to the model with h = 0.1. As observed in pre- 
vious work, there is indeed very good agreement 
between the numerical results for the two models, 
but such an agreement disappears for ^{L)/L > 1, 
where small discrepancies are observed. As it 
can be seen from Fig. 0, the icosahedral points 
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Figure 3. FSS function for the spin-3 susceptibil- 
ity X3- 



Figure 4. FSS function for the spin-4 susceptibil- 
ity X4- 



are above the 0(3) curve and, more importantly, 
the discrepancy increases with L: the points with 
L ~ 10-20 are systematically below the points 
with larger values of L. 

Larger deviations are observed in Figs. || and 
^ where we report the FSS curve for the spin-3 
and spin-4 susceptibilities. Again, the numerical 
results for the icosahedral and the 0(3) model 
agree very nicely up to ^(L)/L ^ 0.8-1, but then 
they indicate that the icosahedral FSS curve is 
steeper than the 0(3) one. Again, notice that 
the discrepancy between the results for two mod- 
els increases with L, indicating that the observed 
effect is not due to corrections to scaling, i.e. it 
is not a lattice artifact disappearing in the con- 
tinuum limit. 

In conclusion, the numerical results we have 
presented show that the icosahedral and the 0(3) 
models belong to different universality classes. 
Note however that discrepancies are observed 
only for ^(L)/i ^ 1, which corresponds to very 
large values of the infinite-volume correlation 
length. For instance, points with L — 128 and 
i{L)lL « 0.94 correspond to ^oo « 10^ 
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